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FOREWORD 


The  work  described  herein  was  conducted  by  the  Mechanics  Section 
of  the  Space  Sciences  Laboratory,  General  Electric  Company,  under  Contract 
NAS  3-8512  with  the  National  Aeronautics  and  Space  Administration,  Lewis 
Research  Center,  Technical  direction  for  the  Lewis  Research  Center  was 
provided  by  S.  Lieblein  and  I.  J.  Loeffler  of  the  Airbreathing  Engines 
Division. 
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ABSTRACT 


The  sensitivity  of  the  peak  axial  pressure  profiles,  obtained  using 
the  impact  model  of  an  earlier  report  (NASA  CR-609),  to  a  change  in  the 
Hugoniot  function  assumed  for  the  reduced  density  projectiles  is  examined 
briefly.  The  previous  Hugoniot  function,  based  on  the  Los  Alamos  equation 
of  state  for  aluminum,  is  replaced  by  one  based  on  the  Plate-Gap  model  of 
a  porous  solid  due  to  Thouvenin.  A  modification  of  this  model  is  described 
in  detail  and  its  comparison  with  existing  experimental  data  for  six  porous 
materials  and  the  predictions  of  the  theoretically  based  equation  of  state  due 
to  Wagner  and  Bjork  for  aluminum  is  noted. 

Comparison  is  also  made  of  calculated  peak  axial  pressures  obtained 
using  the  impact  model  for  the  case  of  normal  density  aluminum  projectiles 
impacting  aluminum  with  experimental  data  obtained  from  corresponding  jet 
projector  impact  shots  conducted  by  the  General  Motors  Defense  Research 
Laboratory. 
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SUMMARY 


The  plate-gap  model  of  a  porous  solid,  due  to  Thouvenin,  is  modified 
by  determining  the  particle  velocity  from  the  momentum-conservation  rela¬ 
tion  and  the  reflected  shock  polar  of  the  non-porous  solid;  the  crush-up  wave 
velocity  formulation  of  Thouvenin  is  retained.  This  procedure  is  shown  to 
yield  Hugoniot  curves  for  solids,  over  a  wide  range  of  initial  porosities, 
which  are  in  good  agreement  with  available  experimental  data  from  Russian 
and  American  sources  for  six  materials. 

The  analytical  model  for  determining  peak  pressure  in  an  impacted 
thick  target,  as  given  in  NASA-CR-609,  yields  a  pressure  profile  in  good 
agreement  with  one  determined  experimentally  by  the  General  Motors  Defense 
Research  Laboratory  for  pressures  above  150  kb  in  aluminum.  Below  this 
level  the  analytical  model,  based  on  hydrodynamic  material  behavior,  is  no 
longer  accurate,  the  deviation  from  the  experimental  profile  becoming  more 
pronounced  as  material  strength  effects  come  into  play. 

Application  of  the  modified  plate-gap  model  Hugoniot  data  to  the 
calculation  of  peak  axial  pressure  resulting  from  impact  of  aluminum  pro¬ 
jectiles  of  varying  degrees  of  porosity  with  a  solid  aluminum  target  resulted 
in  the  removal  of  an  apparent  anomaly  in  the  target  peak  pressure  profile 
obtained  in  a  previous  study  (NASA  CR-609). 
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INTRODUCTION 


In  an  earlier  report  (ref.  1)  an  analytical  formulation  was  presented 
for  evaluating  peak  pressures  in  thick  targets  impacted  by  hypervelocity 
projectiles  of  subnormal  bulk  densities  at  normal  incidence.  It  was  found 
that  the  calculated  initial  peak  pressures  at  the  target  surface  at  impact  were 
greater  for  the  higher  density  projectiles  for  equal  projectile  mass  and 
velocity.  Moreover,  these  pressures  were  never  exceeded  by  those  for  the 
lower  density  equi-energy  projectiles  except  for  the  case  of  the  very  low 
density  projectiles  (0.  1  gm/cc  and  0.2  gm/cc).  For  these  latter  projectiles, 
the  peak  pressure  at  a  given  distance  below  the  target  surface  was  distinctly 
higher  than  that  for  the  projectiles  of  density  0.44  gm/cc  and  above,  (see 
Figures  16,  17  of  ref.  1). 

Since  the  analytical  model  used  in  making  the  pressure  profile  cal¬ 
culations  was  strongly  dependent  on  the  form  of  the  Hugoniot  function  used 
for  the  subnormal  density  projectile  material,  this  anomalous  result  could 
only  be  considered  tentative  until  additional  calculations  could  be  made  using 
other,  perhaps  more  accurate,  Hugoniots  for  describing  porous  material 
under  shock  loading.  The  specific  Hugoniot  used  in  reference  1  was  obtained 
from  the  Los  Alamos  equation  of  state  t  by  eliminating  the  specific 
internal  energy  term  through  use  of  the  energy  jump  across  a  shock  propagat¬ 
ing  in  the  porous  material.  This  Hugoniot,  as  was  pointed  out  in  reference  2 
in  the  case  of  aluminum,  is  not  sufficiently  accurate  for  the  more  highly 
porous  initial  states.  Thus,  it  is  quite  appropriate  to  recalculate  the  pres¬ 
sure  profiles  with  an  alternate  description  of  the  Hugoniot  states  of  a  porous 
solid. 


This  description  is  furnished  by  the  recently  developed  Plate-Gap 
Model  for  a  porous  solid  due  to  Thouvenin  (refs.  3,4).  Thouvenin  replaces  a 
porous  solid  by  a  configuration  of  parallel  plates  of  normal  density,  the 
thickness  and  separation  distance  of  the  plates  being  determined  from  the 
known  degree  of  porosity.  For  this  basic  configuration  Thouvenin  derives 
two  fundamental  relations:  one  for  the  crush-up  wave  velocity  and  the  other 
for  the  particle  velocity  in  the  porous  material  in  terms  of  the  shock  and 
particle  velocities  in  the  normal  density  material.  In  the  present  report, 
Thouvenin's  relation  for  the  particle  velocity  behind  the  crush-up  front  is 
replaced  by  one  which  combines  momentum  conservation  across  the  front  and 
the  assumption  that  the  shocked  porous  material  unloads  along  the  isentropes 
of  the  solid  material,  the  latter  being  taken  as  mirror  images  of  the  Hugoniot 
of  the  solid  in  the  pressure-particle  velocity  plane.  The  plate-gap  model  so 
modified  has  been  used  as  the  basis  for  an  alternate  description  of  Hugoniot 
states  of  an  initially  porous  solid. 

t  This  equation  of  state  is  the  work  of  R.  K.  Osborne  and  his  associates 
in  Group  W-4  at  the  Los  Alamos  Scientific  Laboratory. 
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The  present  report  contains  three  parts.  The  first  part  presents 
the  Plate -Gap  model  in  detail,  including  the  modification  employed  for  cal¬ 
culating  the  particle  velocity,  and  detailed  calculations  of  corresponding 
Hugoniot  data  for  a  wide  range  of  porosities.  Comparisons  are  then  made 
with  available  experimental  data  from  Russian  and  American  sources.  In¬ 
cluded  also  are  Hugoniot  data  for  porous  aluminum  given  by  the  Wagner - 
Bjork  equation  of  state.  The  second  part  compares  calculated  peak  axial 
pressures  obtained  using  the  impact  model  of  reference  1  for  the  case  of 
normal  density  aluminum  projectiles  impacting  aluminum  with  unpublished 
experimental  results  obtained  from  corresponding  jet  projector  impact  shots 
at  the  General  Motors  Defense  Research  Laboratory  (ref.  5).  The  conclu¬ 
ding  section  presents  the  peak  axial  pressure  profiles  for  the  impact  cases 
of  reference  1,  recalculated  using  the  plate-gap  model  Hugoniot  curves,  with 
a  discussion  of  the  differences  which  arise  from  use  of  another  Hugoniot 
representation. 
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THE  PLATE-GAP  MODEL 


In  reference  4  Thouvenin  describes  a  model  for  the  behavior  of  a 
porous  solid  under  shock  loading  which  is  quite  simple.  He  replaces  the 
solid  by  an  array  of  parallel  plates,  each  of  normal  density  pQ  and  of 
thickness  \  pQI  where  PQ  is  the  initial  density  of  the  porous  material; 

the  air  gap  between  any  two  neighboring  plates  is  taken  to  be  1-X  .  (See 

The  porous  material  is  thus  taken  to  be  periodic  in  structure,  the 
width  of  each  period  being  unity  and  of  density  pQ\  =  pQ  ,  thereby  agreeing 
with  the  required  initial  density.  (Symbols  are  defined  in  Appendix  A). 

Consider  now  a  plate  impacting  the  first  plate  at  speed  2u.  A  shock 
is  driven  into  it  at  speed  D(u)  with  corresponding  particle  velocity  u,  where 
D  =  D(u)  is  the  wave  speed  in  the  non-porous  material.  When  the  wave 
reaches  the  rear  free  surface  of  the  first  plate,  the  pressurized  material 
of  the  plate  begins  to  be  unloaded  by  a  rarefaction  wave  moving  back  toward 
the  front  of  the  plate.  The  unloaded  material  then  fills  in  the  gap  1-X  at  the 
free  surface  velocity  ufs  =  2u.  The  total  time  for  all  this  to  happen  is  then 
(see  Figs.  2,3) 
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and  this  is  taken  to  be  the  time  for  the  crush-up  wave  velocity  D*  in  the 
porous  material  to  cover  unit  thickness.  Hence,  we  have  the  first  of 
Thouvenin1  s  relations 


"d*”  =  "d  +  Tu  *  (1) 

wherein  is  used  the  approximation  u  &  2u,  the  validity  of  which  is  discussed 
in  some  detail  in  reference  6.  S 

Thouvenin  derived  a  second  relation  relating  u  and  the  crush-up 

particle  velocity  u*  based  on  an  assumed  equilibrium  crush-up  state.  This 

state  is  assumed  to  occur  when  the  shock  driven  into  the  unloaded  material 

after  its  impact  with  the  second  plate  reaches  an  asymptotic  limiting  strength 

after  being  attenuated  by  the  unloading  waves  moving  to  the  left  in  both  first 

and  second  plates.  Thus,  in  the  time-distance  diagram  of  Figure  3,  the 

shaded  areas  represent  two  neighboring  plates,  the  trajectories  of  the  shocks 

propagating  in  them  at  speed  D  being  given  by  AO  and  BS.  The  trajectory 

of  the  unloaded  material  from  the  rear  of  the  first  plate,  which  fills  the  gap 

1-X  at  speed  2u,  is  the  line  segment  OB  of  slope-2u  with  respect  to  the 

t-axis.  The  crush-up  shock  wave  of  speed  D*  is  represented  by  the  line 

segment  OS  with  slope-D*,  the  ordinate  of  the  point  S  being  t,  =  .  The 

1  D# 


unloading  wave  in  the  first  plate  at  time  t  since  its  initiation  has  the 
trajectory  OM  of  slope  c  Cj  and  Uj  being  the  sound  and  particle 

velocities  in  this  plate  at  this  time.  Similarly,  SM  is  the  unloading  wave 
trajectory  up  to  time  t  from  the  rear  of  the  second  plate  whose  slope  is 
c2-U2  =  (x+1) /(t-tj^).  The  trajectory  of  the  shock  reflected  into  the  unloaded 
material  of  the  first  plate  is  shown  as  the  curve  BM,  whose  slope  at  time  t 
is  the  wave  speed  TJ  =  dx/dt.  Since  a  shock  wave  travels  supersonically 
relative  to  the  material  in  front  of  it  and  subsonically  relative  to  the  material 
behind  it,  the  shock  will  overtake  the  first  unloading  wave  and  will  in  turn  be 
overtaken  by  the  second.  Initially,  U  is  constant  until  the  first  unloading 
wave  is  overtaken,  after  which  U  is  diminished  in  strength  and  the  shock 
trajectory  becomes  curved;  this  attenuation  and  subsequent  curving  is  rein¬ 
forced  upon  the  arrival  of  the  second  unloading  wave  at  time  t  (point  M).  At 
this  time,  Thouvenin  determines  the  trajectory  BM  approximately,  to  within 
second  order  precision,  by  taking  U  to  be  the  average  of  c^-u^  and  c2”u2* 
This  leads  to  the  differential  equation  for  the  shock  trajectory. 


dx  _  _1 
dt  =  2. 


x+1 


t-t 


1  J 


(2) 


whose  solution  is  readily  found  to  be  the  hyperbola 


(3) 


where  k  is  a  constant  of  integration  which  is  determined  by  passing  the  curve 
through  the  point  B  -1,  .  After  some  effort  we  find 


k  = 


2u  +  D*  | 


(l-X  P 
2uX 


(4) 


The  assumed  equilibrium  crush-up  state  is  now  defined  as  that  state  for 
which  U  reaches  its  asymptotic  value,  namely  the  slope  of  the  appropriate 
asymptote  of  (3).  At  this  time  the  speeds  c^-u^  and  c2~u2  will  have  merged 
into  coincidence  with  the  asymptotic  U  value,  so  that 

U  h  c*  -  u*  =  k  -  — - —  ,  (5) 

asymp.  t^ 

where  now  u*  is  the  equilibrium  state  crush-up  particle  velocity  and  c*  is 
the  sound  speed  in  the  material  behind  the  shock  with  strength  given  by  (5). 
Thouvenin  adopts  the  empirical  rule. 
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c  =  A  +  (2B-1)  u  , 


(6) 


relating  sound  and  particle  speeds  behind  a  shock  propagating  in  the  non- 
porous  plate  material,  where  A  and  B  are  empirical  constants,  and  assumes 
the  isentropic  unloading  path  for  the  shock-loaded  porous  solid  (arc  PAB  of 
Figure  2)  to  be  the  mirror  image  of  the  Hugoniot  OP  with  respect  to  the 
vertical  PC.  With  this  assumption  c*  can  be  obtained  by  replacing  u  in  (6) 
by  2u  -  u*,  whence  we  find 


c*  -  u*  =  A  +  (2B-1)  2u  -  2B  u*  . 
Combining  (5)  and  (7)  then  yields  the  following  result  for  u*: 


u* 


1 

2B 


A  +  2  (2B-1)  u+D*  -k  , 


(7) 


(8) 


where  u  is  known  and  D*  and  k  are  obtained  from  (1)  and  (4),  respectively. 

The  constants  A,  B  are  the  same  ones  which  appear  in  the  classical  empirical 
relation,  D  =  A  +  Bu,  relating  shock  and  particle  velocities  in  a  given  ma¬ 
terial,  for  which  an  extensive  tabulation  is  given  in  reference  6. 


Equations  (3),  (4),  (5),  (8)  do  not  appear  in  Thouvenin's  paper  but  are 
derived  here  on  the  basis  of  the  procedure  which  Thouvenin  prescribes. 
Actually,  Thouvenin  simplified  the  procedure  still  more  by  omitting  one  of 
the  unloading  waves  in  Figure  3,  thereby  simplifying  the  differential  equation 
(2)  and  obtaining  a  parabola  for  the  shock  trajectory.  He  does  not  give  the 
corresponding  formula  for  u^,  However,  he  does  plot  p^  vs  u for  porous 
copper  at  moderate  pressures,  where  relation  (8)  agrees  with  his  plot 
reasonably  well.  At  higher  pressures,  as  might  be  expected  from  the  type 
of  approximations  made,  relation  (8)  is  quite  inaccurate.  For  this  reason 
an  alternate  relation  for  the  computation  of  u*  to  replace  equation  (8)  has  been 
developed. 
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MODIFIED  PLATE-GAP  MODEL 


The  experimental  evidence  presented  by  Thouvenin  (ref.  4)  in  support 
of  relation  (1)  for  the  case  of  porous  copper  and,  to  a  lesser  extent,  porous 
aluminum  and  uranium.,  is  most  convincing.  For  this  reason  relation  (1)  will 
be  assumed  to  be  valid  without  amendment.  It  is  proposed,  however,  to 
calculate  u*  by  the  following  procedure.  One  notes  from  the  physics  of  the 
situation  and  from  Figure  2  that  u*  will  lie  between  u  and  2u  and  that  by 
virtue  of  symmetry  the  corresponding  pressure  p*  can  be  obtained  from  the 
shock  polar  for  the  non-porous  solid,  p  =  F(u),  evaluated  at  2u-U'‘'.  Thus, 

p*  =  F  (2u-u*)  .  (9) 

If  now  one  keeps  the  Hugoniot  relation  for  conservation  of  momentum  across 
the  crush— up  shock. front,  we  shall  also  have 

p*  =  Po  D*  u*  .  (10) 

From  the  equality  of  (9)  and  (10),  we  then  have  the  necessary  relation  for 
calculating  u*: 

F  (2u-u*)  =  p  D*u*  (11) 


This  is  in  contrast  to  Thouvenin’s  procedure  of  calculating  u*  by  an  inde- 
pendent  relation  like  (8)  and  then  getting  p*  from  (9). 

For  a  given  material,  then,  of  porosity  X  the  modified  procedure  for 
calculating  shock  velocity  D*,  induced  by  plate  impact  at  speed  2u,  and  the 
pressure  p*  behind  this  shock  with  accompanying  particle  speed  u*,  is  sum¬ 
marized,  respectively,  by  relations  (1),  (10)  and  (11). 

Calculated  curves  of  p*  and  D*  vs  u*,  and  of  p*  vs  p*/PQ  for 
porosities  m  ranging  from  1  to  6.  14  are  presented  in  Figures  4  through  9 
for  aluminum,  copper,  lead,  iron,  nickel  and  tungsten.  The  (D--',  u-')-plots, 
for  given  X>  are  based  on  equations  (1),  (11)  with  u  regarded  as  a  parameter, 
the  function  F(u)  for  the  non-porous  solid  was  available  from  the  Los  Alamos 
formulation  (see  pp.  13,  14  of  ref.  1).  Values  of  D*  and  u*.  when  substituted 
into  equation  (10),  furnish  values  of  p*  for  the  (p*,u*) -plots.  The  values  of 
p*/p  to  be  coupled  with  the  p*  values  are  given  by  the  continuity  of  mass 
relation  across  the  crush-up  front,  namely 


XD* 

D--u*  * 


(12) 
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Experimental  data  for  each  of  the  six  metals,  taken  from  references  6 
through  11,  are  also  plotted.  In  addition,  in  the  case  of  aluminum,  points 
(Fig.  4)  are  plotted  from  the  Hugoniot  curves  implied  by  the  Wagner-Bjork 
equation  of  state  (ref.  2). 

As  is  evident,  the  modified  plate-gap  model  computations  show 
excellent  agreement  with  experiment  and,  in  the  case  of  aluminum,  with  the 
Wagner-Bjork  points,  on  the  (p*.  u*)  and  (D*,  u*)  plots.  On  the  plots  of 
pressure  versus  compression,  (p*,  p*/pQ),  the  agreement  between  theory 
and  experiment  ismixed:  excellent  for  some  of  the  data  points,  poor  for 
others.  In  view  of  the  excellent  agreement  in  the  case  of  the  (D*,  u*)  and 
(p*,  u:)- plots,  the  reason  for  this  is  not  altogether  clear.  Experimentally, 
it  is  D*  that  is  measured  and  the  other  variables  are  then  calculated  from  it 
using  the  Hugoniot  relations.  The  relative  errors  in  p*  and  u*  are  related 
to  the  relative  error  in  D*  by  (see  Appendix  B) 


A  u*  _  _1  AP*  Ap*  1  A  D* 

u*  2  D*  *  p*  “  2  D*  *  (1^) 

Thus,  a  good  measurement  of  D*  yields  smaller  relative  errors  for  p*  and 
u*.  The  relative  error  in  the  compression,  cr*  =  p*/p  ,  as  equation  (B.  6) 
shows,  depends  upon  the  compression  and  the  porosity:  ° 


a* 


3 

2 


(ma*  -  1) 


A  D* 
D* 


(14) 


2 

The  factor  -  (ma*-l)  is  larger  than  1  whenever  a*  >  —  .  The  latter  was 

true  for  all  but  one  of  the  experimental  porous  data  points  plotted.  This 
suggests  that,  in  general,  relative  errors  in  the  compression  are  larger 
than  corresponding  relative  errors  in  the  measured  shock  velocity  D*,  the 
amplification  increasing  with  the  porosity.  The  departure  of  the  experimental 
data  points  from  the  modified  plate -gap  model  curves  in  the  pres  sure -com¬ 
pression  plane  may  be  due  both  to  the  errors  in  the  pressure  and  the  com¬ 
pression  resulting  from  a  given  error  in  the  shock  velocity  and  also  to.  the 
fundamental  inadequacy  of  the  model  itself.  How  to  separate  these  effects 
m  evaluating  the  model  is  an  additional  problem  for  which  no  answer  is  now 


APPLICATION  TO  IMPACT  BY  REDUCED  DENSITY  PROJECTILES 


An  analytical  model  for  determining  peak  axial  pressure  in  a  thick 
target  impacted  at  normal  incidence  by  end -on  oriented  right  circular  cylin¬ 
drical  projectiles  of  either  normal  or  subnormal  bulk  density  travelling  at 
hypervelocity  was  presented  in  detail  in  reference  1.  This  model  takes  into 
account  the  speed  and  geometry  (radius  L  and  length  i)  of  the  projectile,  its 
density  and  that  of  the  target,  and  the  speed  of  sound  in  the  undisturbed  tar¬ 
get  material.  It  is  based  on  combining  the  Hugoniot  relations  for  the  ther¬ 
modynamic  state  back  of  the  shock  advancing  into  the  target  with  a  kinematic 
relation  expressing  the  axial  speed  of  the  shock  as  a  function  of  its  axial 
position.  Empirical  constants  appearing  in  the  kinematic  relation  are  evalu¬ 
ated  from  a  knowledge  of  the  location  of  the  axial  point  where  the  shock  is 
first  attenuated  by  lateral  rarefaction  waves  originating  at  the  projectile 
periphery  at  impact  and  from  empirical  formulations  based  on  numerical 
solutions  of  the  equations  defining  the  impact  generated  flow  in  the  target. 

The  kinematic  relation  itself  is  based  on  a  blast  wave  interpretation  for  the 
axial  shock  speed,  expressing  the  latter  as  a  linear  combination  of  the  speeds 
of  planar  and  cylindrical  blast  waves. 

Experimental  data  are  now  available  to  verify  the  accuracy  of  the 
model  (ref.  5).  These  data  of  pressure  against  distance  in  aluminum  tar¬ 
gets  were  obtained  by  the  General  Motors  Corporation  Defense  Research 
Laboratory  under  contract  with  the  NASA-Lewis  Research  Center.  The 
experimental  data  are  presented  for  an  impact  into  an  1100-0  aluminum 
target  at  11.  15  km/ sec.  by  a  solid  aluminum  projectile  weighing  0.  585  gm 
with  a  length  to  diameter  ratio  (X  /2L)  of  3.  3.  In  reference  5  peak  axial 
pressure  against  distance  into  the  target  was  represented  by  a  straight 
line  on  log-log  paper.  This  relation  is  shown  by  the  solid  line  in  Fig¬ 
ure  10. 


When  the  peak  axial  pressure  model  of  reference  1  is  applied  to  this 
particular  impact  case,  the  profile  obtained  is  the  curved  dashed  line  of 
Figure  10.  Above  150  kb  this  profile  is,  within  experimental  error,  coin¬ 
cident  with  the  experimental  profile  of  reference  5.  Below  150  kb  the  two 
profiles  diverge  since  the  profile  of  reference  1  becomes  invalid  at  lower 
pressures.  This  is  because  the  profile  of  reference  1  is  based  on  a  purely 
hydrodynamic  model  of  material  behavior.  Below  150  kb  the  aluminum  tar¬ 
get  begins  to  exhibit  behavior  markedly  different  from  that  of  a  perfect 
fluid. 

In  reference  1,  hypervelocity  impacts  of  equi-mass  aluminum  pro¬ 
jectiles  of  densities  2.  702,  0.9,  0.44,  0.2  and  0.  1  gm/cc  into  a  thick 
aluminum  target  were  investigated  using  the  peak  pressure  model.  Two 
impact  velocities  were  considered,  20  km/sec  and  7.6  km/sec.  The  purpose 
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was  to  investigate  variation  in  peak  axial  pressure  propagating  into  the  tar¬ 
get  as  a  function  of  projectile  density.  It  was  found  that  at  20  km/sec,  the 
peak  pressure  profiles  in  the  target  due  to  impacts  of  the  projectiles  of  den¬ 
sities  2.  702,  0.  9  and  0.44  gm/cc  were,  beyond  a  short  distance  into  the 
target,  essentially  indistinguishable.  For  the  low  density  projectiles  of  0.  1 
and  0.2  gm/cc,  however,  the  peak  pressure  profiles  generated  in  the  target, 
beyond  a  certain  distance  into  the  target,  lay  above  the  profiles  arising  from 
the  higher  density  projectile  impacts.  The  same  situation  was  found  to  hold 
for  the  profiles  generated  by  impacts  at  7.  6  km/sec.  As  was  pointed  out  in 
reference  1,  this  anomalous  result  could  only  be  considered  tentative  in 
view  of  the  strong  dependence  in  the  analytical  model  on  the  form  of  the 
Hugoniot  function  used  for  the  subnormal  density  projectile  material.  The 

situation  will  now  be  re-examined  using  an  alternate  Hugoniot  representa¬ 
tion. 

For  the  peak  axial  pressure  model  of  reference  1,  the  Hugoniot 
values  of  pressure,  density,  shock  and  particle  velocities  in  the  target  at 
impact  are  required  as  inputs.  In  Figures  16  and  17  of  reference  1,  the 
Hugoniot  values  were  obtained  based  on  the  Los  Alamos  equation  of  state  for 
aluminum.  The  present  calculations  use  Hugoniot  values  based  on  the 
Modified  Plate -Gap  Model. 

Figures  11  and  12  show  the  plots  of  recalculated  peak  axial  pressure 
versus  shock  penetration  distance  into  a  solid  aluminum  target  arising  from 
impacts  at  speeds  of  20  and  7.6  km/sec,  respectively,  of  four  aluminum 
projectiles,  all  of  the  same  mass,  but  of  densities  2.702,  0.90,  0.44,  and 
0.20  gm/cm  .  The  projectiles  are  taken  to  be  right  circular  cylinders  with 
length  l  and  radius  L  equal  to  1/2  ,6  and  are  assumed  to  impact  end -on  at 
normal  incidence.  The  values  of  L  associated  with  these  densities  are, 
respectively,  0.2619,  0.  3778,  0.4797,  and  0.6238  cm. 

Figure  11  is  noticeably  different  than  Figure  16  of  reference  1. 
Whereas,  formerly,  the  pressure  profiles  for  reduced  density  projectiles 
of  0.44  gm/cm3  and  higher  showed  coalescence,  indicating,  at  least  for 
this  range  of  densities,  no  pronounced  density  effect,  the  profiles  corres¬ 
ponding  to  the  reduced  density  projectiles  now  lie  distinctly  below  that  for 
the  projectile  of  normal  density.  In  addition,  the  formerly  anomalous 
result  in  whij:h  the  profile  corresponding  to  the  projectile  of  density 
0.  20  gm/cm  lay  above  the  normal  density  profile  is  no  longer  present. 

If  one  interprets  cratering  in  a  thick  target  in  terms  of  the  impact¬ 
generated  pressure  pulse  propagating  therein,  the  anomaly  of  reference  1 
suggests,  contrary  to  physical  intuition,  that  very  low  density  projectiles 
have  greater  penetration  power  into  thick  targets  than  higher  density  pro¬ 
jectiles  of  the  same  mass  and  impact  velocity.  The  present  result,  for 
20  km/sec  impacts,  implies  then  reduced  crater  depths  in  thick  targets 
impacted  by  highly  porous  projectiles. 


10 


An  apparent  anomaly  is  still  present  in  Figure  12  for  the  7.  6  km/ sec 
impact.  The  pressures  in  Figure  12  for  the  0.20  gm/cm3  case  all  lie  below 
100  kb,  and  in  this  low  range  of  pressures  the  peak  axial  pressure  model  of 
reference  1  no  longer  applies  since  it  was  based  on  a  purely  hydrodynamic 
model  of  material  behavior.  Credence  for  this  statement  is  supplied  by 
Figure  10,  where  an  experimental  pressure  profile  for  an  impact  into 
aluminum  of  a  solid  aluminum  projectile  of  rra  ss  0.  585  gm  at  11.15  km/ sec 
is  compared  with  that  predicted  by  the  model  of  reference  1.  The  two  pro¬ 
files  coincide  from  1  Mb  down  to  150  kb,  after  which  there  is  rapid  diverg¬ 
ence.  Thus,  it  is  argued  that  for  the  7.  6  km/ sec  impact  of  a  projectile  of 
density  0.20  gm/cm3,  the  location  of  its  pressure  profile  above  that  of  the 
normal  density  profile  is  meaningless  since  at  the  low  pressures  involved 
neither  of  the  profiles  is  valid. 

The  conclusion  which  appears  warranted  from  these  results 
is  that  the  anomalous  inversion  of  the  profiles  corresponding  to  projectile 
densities  of  0.20  and  2.  702  gm/cm3  for  equi-energy  impacts  at  20  km/sec 
arises  from  the  use  of  inaccurate  Hugoniot  data  for  the  reduced  density 
projectile  and  that  the  same  situation  at  7.  6  km/sec  is  invalid  owing  to  the 
inapplicability  of  the  analytical  model  at  the  low  pressures  encountered. 

As  demonstrated  by  the  curves  in  Figures  4  thru  9  with  relation  to  the 
available  experimental  data,  the  modified  plate -gap  model  can  be  relied 
upon  to  furnish  accurate  Hugoniot  data  for  impact  calculations. 


CONCLUSIONS 


1.  The  plate-gap  model  of  a  porous  solid,  as  modified  in  this  report, 

can  be  used  to  generate  Hugoniot  data  for  solids  over  a  wide  range  of  initial 
porosities.  This  is  confirmed  by  good  agreement  with  available  experimental 
data  from  both  Russian  and  American  sources, 

^  *s  shown  that  the  analytical  model  for  determining  peak  axial 
pressure  in  an  impacted  thick  target,  as  developed  in  NASA  CR-609,  yields 
a  pressure  profile  in  good  agreement  with  one  determined  experimentally 
by  the  General  Motors  Defense  Research  Laboratory  for  pressures  above 
150  kb  in  aluminum.  Below  this  level  the  analytical  model,  based  on  hydro- 
dynamic  material  behavior  is  inapplicable. 

3.  When  the  modified  plate-gap  model  Hugoniot  values  are  applied  to  the 

calculation  of  peak  axial  pressures  resulting  from  impact  of  aluminum  pro¬ 
jectiles  of  varying  degrees  of  porosity  with  a  solid  aluminum  target,  pre¬ 
viously  obtained  results  of  an  anomalous  nature  (NASA  CR-609)  either 
disappear  or  can  be  shown  to  be  inapplicable.  At  a  given  axial  position  in 
the  target,  reduced  projectile  density  results  in  reduced  pressure  in  all 
applicable  cases. 
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APPENDIX  A 


SYMBOLS 

A  constant  in  relation  D  =  A  +  Bn 

B  constant  in  relation  D  =  A  +  Bn 

c  speed  of  sonnd 

D  shock  velocity 

F  shock  polar  function 

L  radins  of  cylindrical  projectile 

l  length  of  cylindrical  projectile 

m  porosity  ratio,  p  /p 

o  o 

p  pressure 

t  time 

n  particle  velocity 

x  abscissa 

X  1/m 

p  density  behind  shock 

p  undisturbed  target  density 

o 

p  reduced  density  of  porous  projectile 

o 

Subscripts 

fs  free  surface  value 

Superscript 

*  indicates  value  of  given  quantity  for  porous  material 


APPDENDIX  B 


DETERMINATION  OF  RELATIVE  ERROR  RELATIONS 


For  a  planar,  steady-state  shock  of  strength  p*  propagating  in  a 
porous  material,  the  work  performed  per  unit  area  per  unit  time  is  constant, 
that  is 


p*  u*  =  const. 


Since  p*  =  Xp  D*  u*, 
o 


we  have 


D''c  u* 


=  const, 


whence 


An*  _  I  AP* 

u*  “‘2D*’  (B«  1) 

From  the  conservation  of  mass  across  the  crush-up  front,  denoting 
P*/pQ  by  cr*,  we  have 


u*  mcr*-l 

D*  mg*  *  (B*2) 


whence 


Au *  =  AD*  1  Aa* 
u*  D*  ma*-l  a  * 


From  p*  =  Xp  D*u*,  we  have 
o 


A  p*  _  A  D*  A  u* 
p*  “  D*  +  u* 


Combining  (B.  1)  and  (B.4)  then  yields 


A  p*  _  _1  A  D* 
p*  2  D* 

Similarly,  combining  (B.  1)  and  (B.  3)  yields 


(B.3) 


(B.4) 


(B.5) 


Act  * 

cr  * 


3 

—  (ma  *  - 1 ) 


AD* 

D*  ' 


(B.  6) 
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Figure  4.  Hugoniots  for  Porous  Aluminum  -  Modified  Plate  - 
(a)  Shock  Velocity  vs  Particle  Velocity 
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Figure  4.  Hugoniots  for  Porous  Aluminum  -  Modified  Plate-Gap  Model 
(b)  Pressure  vs  Particle  Velocity 
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Figure  5.  Hugoniots  for  Porous  Copper  -  Modified  Plate-Gap  Model 
(a)  Shock  Velocity  vs  Particle  Velocity 


PRESSURE,  p*(mb) 


=  6.14 


Hugoniots  for  Porous  Copper  -  Modifi 
(c)  Pressure  vs  Compression 


SHOCK  VELOCITY,  D7(CM//j.SEC) 


Shock  Velocity  vs  Particle  Velocity 
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Figure  8.  Hugoniots  for  Porous  Nickel  -  Modified  Plate-Gap  Model 
(b)  Pressure  vs  Particle  Velocity 
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Figure  9.  Hugoniots  for  Porous  Tungsten  -  Modified  Plate-Gap  Model 
(a)  Shock  Velocity  vs  Particle  Velocity 
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Figure  9.  Hugoniots  for  Porous  Tungsten  -  Modified  Plate-Gap  Model 
(b)  Pressure  vs  Particle  Velocity 


I  I  I  I  I 


o 

—  oc 


— :  I — 
LiJ 
O 
Q C 
00  <t 


< 


PO  CVI I  “  O  00.  (£>  lf>  ;  IO  cvl 


^  <p 

*  s®  • 

rt  ^  u 

(D  <w  (p 
p_,  H  OT 

xf  g  "d 
$  §  & 
.—4  *d  VO 

2  S  • 

V  3  t>- 

3«  * 


(qi^)*d  ‘3ynSS38d  1VIXV  XV3d 


NASA-Langley,  1968  - 


